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Abstract. - Quantum trajectories defined in the de Broglie-Bohm theory provide a causal 
way to interpret physical phenomena. In this Letter, we use this formalism to analyze the short 
time dynamics induced by unstable periodic orbits in a classically chaotic system, a situation 
in which scars are known to play a very important role. We find that the topologies of the 
quantum orbits are much more complicated than that of the scarring and associated periodic 
orbits, since the former have quantum interference built in. Thus scar wave functions are 
necessary to analyze the corresponding dynamics. Moreover, these topologies imply different 
return routes to the vicinity of the initial positions, and this reflects in the existence of different 
contributions in each peak of the survival probability function. 



Trying to overcome interpretative difficulties T encountered in the standard quantum 
theory Bohm developed in the 1950's an alternative formalism, based on the notion of 
quantum trajectories [3]. This complementary quantum theory of motion ^ combines the 
accuracy of the standard quantum description with the intuitive insight derived from the 
causal trajectory formalism, thus providing a powerful tool to understand the physical mech- 
anisms underlying microscopic phenomena. A key property of this theory is non-locality, 
since particles are guided along their paths by a pilot wave, constituted by the system wave 
function 0. 

The above mentioned characteristics are the cause of the vigorous revival that the method 
is presently experiencing. Taken from its computational side, it has been succesfully applied 
to treat different realistic physical systems 6 , 810111. This way of proceeding allowed, for 
example, Wyatt ^U] to analyze tunneling of a wave packet through a barrier, nicely elucidating 
the corresponding mechanism in terms of quantum trajectories. A portion of them, in the 
foremost part of the initial packet, overcome the obstacle because they acquire some extra 
kinetic energy during a boost phase, shortly after the launching of the particle. With a 
similar analysis in mixed quantum-classical simulations, Prezhdo and Brooksby solved 
the quantum backreaction problem for the dissociation reaction of O2 on Pt. 
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Since the de Broglie-Bohm (BB) theory considers trajectories in a fuUy quantum framewok, 
another hne of research to which it has been naturaUy apphed is quantum chaos |12ll3j . Using 
this idea, several authors '14' have foUowed the separation of nearby Bohmian trajectories, thus 
computing quantum analoges to Lyapunov exponents. In our opinion, this way of proceeding 
is too simple. The relative behavior of two particular trajectories do not mean, in general, 
anything in view of the hoUistic character 4 imposed by the non-local character of BB 
theory 9^. Moreover, the quantum Lyapunov exponents so defined are quantities introduced 
in the context of infinitely long times, while it is widely recognized that many interesting 
quantum phenomena, such as scars 15 , have its origin in the short time dynamical evolution 
of wave packets, much before the so called log time is reached. The relevant information 
about short time recurrences, and other quantum time scales, is contained in the survival 
probability, ^(t) = 1(0(010(0)) |2. 

In this Letter we present a study of the topology and characteristics of quantum trajectories 
computed using the BB theory in a situation in which the dynamics of the initial state is 
strongly determined by a classically chaotic PO. The idea is to take advantage of the intuitive 
character of the BB mechanics in order to provide a good understanding of the quantum short 
time dynamics taking place in the process, complementary to that provided by the usual scar 
theory [T2[TS1|THI. 

The fundamental equation in the BB theory is derived from Madelung hydrodynamical 
formulation of quantum mechanics 17 , which recasts the wave function in polar form, 

0(r,i) = i?(r,t) e'^^""^*) (1) 

(throughout the paper h is set equal to 1). Quantum trajectories are then defined by means 
of expression 

mf ^\/S^ Im(V0/0), (2) 

which can be numerically integrated once (f> is known. Alternatively, one can consider the cor- 
responding newtonian form, mr = —V{V + Q), where the extra term, Q = — (l/2m)(V^i?/i?) 
called quantum potential, is responsible for all quantum effects taking place in the process. 

The problem that we have chosen to study is the dynamics of a particle of mass 1/2 
enclosed in a desymmetrized stadium billiard of radius r = 1 and area 1 + 7r/4 with Dirichlet 
boundary conditions, influenced by the diagonal periodic orbit (PO). This corresponds to a 
diamond shape in the full version of the stadium. For this purpose, we start from a harmonic 
oscillator coherent state 

0(0) = (2aA)i/4e-"(— o)^-a(y-yo)^+i(P^+P„"y)^ (3) 

with a = 30.68, initially located on the middle of the diagonal, {xq, yo,P^,Py) = (1,1/2, 96/-\/5, — 48/a/5). 
From these values we have E — 2304 and a period T = 0.0466. As stated in the introduc- 
tion, the associated short time dynamics can be adequately discussed in terms of the survival 
probability, S{t), which is shown in full line in Fig.QJup to t = 0.1. As can be seen, the curve 
consists of two peaks, taking place at t ~ 0.047 and ^0.094, respectively, originated by the 
first two recurrences of the initial packet ©. Also, the second peak shows some substructure, 
in the form of two shoulders located to the left of the main maximum. These features were 
explained jJSj in terms of the scarring effect of the dynamics along the diagonal PO, using 
dynamically averaged wave packets 19;. With this method, we were able to uncover the role 
played by other short, but longer, POs dynamically related to the diagonal one. In the inset 
to Fig. ^we have represented the diagonal PO, contributing to both recurrences in S{t), plus 
the extra PO contributing also to the second one. 
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Fig. 1 - Survival probability function corresponding to a wave packet initially centered on the diagonal 
periodic orbit of a desymmetrized stadium billiard with r = 1 and area l+7r/4, calculated both exactly 
(full line), and with the quantum trajectories approximation defined in the text (dashed line). The 
insets show the two periodic orbits responsible for the two recurrence peaks observed. 



In this Letter we want to go step further and provide new insight in these results, using 
the interpretative abilities of the BB theory to understand quantum effects ^J^j. For 
this purpose, we will use an approximate version of the survival probability, defined in the 
following way. Starting from the phase space representation in terms of the Wigner function, 
W, associated to 0, 



S{t) 



J J dqdPM^o(q,P)m(q,P), 



(4) 



W is substituted by the following approximate distribution based on quantum trajectories, 
[Xi{t),y,{t)], 



Pt{x,y,P.^,Py)=N 



N 



5[x - Xi{t)] 5[y - yi{t)] 



x6 



P. 



dS_ 

dx 



p, 



dy 



y=yi{t). 



(5) 



Finally, and in order to make p a smooth function, gaussian distributions on the delta functions 
are allowed, thus obtaining 



N 



Sit) « iV-1 J2 exp{-a[x,(t) - x,(0)]2 - a[y,it) - y,{0)f 



(6) 



The results obtained with this approximation are shown in dashed line in Fig. ^ They 
have been computed using 80 quantum orbits, obtained by numerical integration of eq. jSJ, 
for which a value of cr = 156.25 is adequate. Some caution must be exerted when integrating 
the quantum equations of motion in cases like ours, for which the dynamics of the system are 
very chaotic. As will be seen below, the time-dependent wave packet develops very quickly 
a tremendous complexity, showing a complicated nodal structure of size of the order of fi. 
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Fig. 2 - Propagation of 80 quantum trajectories corresponding to a wave packet initially centered on 
the diagonal periodic orbit of a desymmetrized stadium billiard with r = 1 and area 1 + n/4. See 
text for details. 

Fig. 3 - Time evolution of the probability density, \(f){x,y,t)\^ , corresponding to the quantum trajec- 
tories of Fig. |5| See text for details. 



In the present case we have chosen a Gear stiff integration method with tolerance control. 
The trajectories were initially distributed uniformly over four radii: 0.00333, 0.0333, 0.05 
and 0.1 respectively, around {xo,yo)- As can be seen in Fig. the agreement between the 
global shape of the exact and approximate survival probabilities is good, except for some 
discrepancies in the relative intensities. Moreover, both calculations present recurrences at 
exactly the same times, and the latter even reproduces the shoulder structure existing in the 
exact second recurrence. It should be stressed at this point that our approximate results can 
be easily converged to the exact ones, by using a better (gaussian) sampling procedure of the 
initial conditions. A full account of this treatment will be given elsewhere [201 • However, this 
is not an essential point for the discussions of this Letter, in which we prefer, for the sake 
of the discussions below, to use our imperfect sampling, that artificially hide the importance 
of quantum trajectories arriving to the first recurrence at average times, in favour of those 
making it either too early or too late. Due to this fact, our approximate results (see Fig. 
clearly show the existence of two contributions, which nicely reproduce both tails of the first 
recurrence. 

The fact that we are basing our analysis in an approximate expression may be considered 
at first sight as a drawback. However, it is quite the contrary, since in the way in which it is 
defined, our version of S{t) comes given as a sum of positive contributions, this having the 
important advantage of allowing the localization of those trajectories contributing the most to 
S{t). This makes our method a superb tool to investigate, in an intuitive and straighforward 
way, how the dynamics of the initial wave packet takes place. 

To discuss this in detail, let us start by examining the quantum trajectories used in our 
calculation of S{t), which are shown in Fig. El In order to make the figure clear, these 
trajectories have been splitted into several consecutive segments, which are presented in the 
different panels. In this way, the first plot, labeled tf = 0.012, corresponds to paths described 
by the trajectories during interval t = — 0.012, the second in interval t = 0.012 — 0.023, and 
so on. To help in the discussion, it is interesting to consider also the time evolution of the 
corresponding wave packet, which is shown in Fig. |3| These results, that were used in the 
integration of eq. El have been obtained by projecting 0(0) on the stadium eigenstates and 
subsequent application of the evolution operator. 
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Fig. 4 - (Upper part) Quantum trajectories representing the largest contributions to the approximate 
version of S{t) [eq. |S| for peaks (a), (b) and (c) in Fig. 

(Lower part): Scar wave function, constructed with the average method of Ref. TST, corresponding 
to trajectories in (a) and (b) (left part), and to (c) (right part). 

As expected, trajectories (Fig.O and wave packet (Fig.|3J) follow, in the first part ot their 
evolution {t < 0.023), a very classical path, as dictated by Ehrenfest's theorem. After the 
first rebound at the lower right corner, they both spread at t ~ 0.025 in a fan-type pattern, 
experiencing the well known focalization effect originated by the self-focal point of the orbit 
|16| . Afterwards, the dispersed wave packet (see Fig. I^J collides with the region around the 
upper left corner, giving rise to a noticeable series of horizontal fringes, formed by the maxima 
and nodal lines of the wave function. These maxima originate fiat plateaus in the associated 
quantum potential, separated by deep canyons along the nodal lines. This pattern greatly 
affect the topology of the corresponding quantum trajectories, as can be seen in Fig.|21 Then, 
trajectories, during the interval t ~ 0.034 — 0.056, are mostly horizontal, following paths along 
the maxima of the corresponding pilot wave (see central tier in Fig. I^J- Moreover, particles 
are "kicked" when trying to cross this structure, giving rise to the short diagonal segments 
observed, for example, in some trajectories at the top of panel labelled t = 0.034. Other 
rebounds take place at subsequent times, originating for t > 0.06 a complicated structure 
in the distribution of the quantum probability density. The topology of the corresponding 
trajectories becomes accordingly also very complex, showing an extreme sensitivity to both 
initial conditions and the accuracy imposed to the integration procedure. It should, however, 
be pointed out the extreme resemblance which is obtained between the patterns formed by 
the maxima of (j){t) in Fig. El and the accumulations of quantum trajectories in Fig. |21 This is 
even more surprising when considering the small number of trajectories used in our calculation. 
The above mentioned dynamical sensitivity can be a problem when trying to interpret the 
behavior of individual trajectories for very long times, but it is reasonable to admit that the 
statistical properties of ensembles of trajectories, such as that considered in the approximate 
calculation of S(t) [eq. ©], behave adequately, similarly to what happens in ergodic theory. 

After having discussed the general behavior of the quantum trajectories used in our study, 
let us investigate now which one of them are more relevant to explain the different features 
exhibited by the survival probability. Figure ^ shows such trajectories for the three labeled 
peaks of S{t) in Fig. Q Panels (a) and (b) corresponds to trajectories in the two separated 
peaks at the first recurrence. Those in (a) correspond to particles that arrive earlier to the 
vicinity of the initial distribution, p(0), thus corresponding to what can be called forerun- 
ners. As can be seen, the topology of the associated trajectories is extremely similar to the 
classical diagonal PO, showing only a very small dispersion, that mirrors the variation of the 
corresponding initial conditions. 
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On the other hand, trajectories in panel (b), arriving later to the initial region, trace 
slightly different paths. At the upper left corner, they make a left turn that delays their arrival 
to the final destination. By doing this, they allow the pilot wave to spread more profusely, 
also giving it enough time to develop the horizontal series of nodal lines that are observed in 
Fig. As explained before, the existence of these nodal lines greatly affects the paths taken 
by quantum trajectories, being the origin of the ladder-type structure observed in the last part 
of trajectories shown in Fig.^b). This important quantum effect has never been discussed in 
the literature before, and shows that the peaks in S{t) consist of different contributions with 
somewhat heterogeneous histories, that determine the time delay with which they arrive to 
the initial region. This is a nice example of how Bohmian mechanics can deal with the time 
variable in a straighforward way |21j . in contrast to what happens in the standard quantum 
theory. Also, these contributions can be traced back to their original launching positions in 
the initial wavepacket, thus establishing a true causal connection between initial and final 
conditions |20| . 

Let us analyze next the second recurrence at i ~ 0.09. Now the situation is much more 
complicated, since as described in Ref. jTH] at least two other POs need to be considered. 
In view of our approximate results, the most interesting feature to considered here is the 
maximum, leaving aside the shoulder structure. Figure ^c) shows the quantum trajectories 
contributing the most to that part of S{t). As can be seen, an attempt to give an interpretation 
in the same terms used before for the first recurrence is less convincing, since the existence 
of a connection between quantum trajectories and lines of maxima at the corresponding time 
(see panels in Fig. |3J) is not so obvious. This is not unexpected since, as stated before, the 
dynamics here get diluted into too many POs. 

An alternative way to interpret these results is to use scar functions. This is going one 
step further the classical dynamics of the POs, which only represents a first approximation to 
the problem. In quantum (or semiclassical) mechanics one has to include (at least to some 
extend) the associated wave properties. Scar wavefunctions can be constructed in a number 
of ways. We will use the method of dynamical averaging introduced by us 19,22 , that is very 
efficient. A very elegant way to construct analogous functions has been discussed by Vergini 
and Carlo in their resonance treatment of quantum chaos |23| . The scar wave function relevant 
for the discussion of the first recurrence in S{t) is presented in the left part of Fig. 01 below 
the corresponding trajectories [panels (a) and (b)]. As can be seen, it presents the right 
physical characteristics to explain the (slightly) different behavior of the two peaks uncovered 
by our approximate calculation of the survival probability. Namely, the function is on the one 
hand localized along the diagonal PO, with a series of nodal lines perpendicular to it, and 
a distribution of the probability density in agreement with the expected focalization effect. 
But more importantly, this function presents at the upper left corner the series of horizontal 
fringes explaining the delayed behavior of the trajectories contributing to the second peak of 
the first recurrence in S{t) [trajectories in (b)]. 

Finally, with respect to the maximum of the second recurrence the corresponding scar wave 
function is shown in right part of Fig.0] In terms of it, the interpretation of trajectories in (c) 
is much more obvious, and analogous to that given before. The relevant quantum trajectories 
clearly follow the lines of maxima of the associated scar wave function. The reason for the 
agreement is clear. The construction procedure of the scar function only includes dynamical 
information concerning an isolated PO and its associated manifolds "5?, thus making possible 
the comparison that we were looking for. 

In summary, we have presented in this Letter a first study on the dynamical behavior of 
quantum trajectories, as defined in the de Broglie-Bohm theory, for a system in which the 
classical dynamics are highly chaotic. Searching for the elusive concept of quantum chaos, this 
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situation had been previously studied by other authors in the infinite time hmit. Contrary, we 
have considered here the other Hmit, i.e. short time dynamics, for which scarring effects are 
known to be important. Proceeding in this way and taking advantage of the causal character of 
our approach we have found new interesting features in the corresponding survival probability. 
In particular, we have shown how the peaks appearing in this function at the recurrence times 
can be decomposed in different contributions with heterogenous dynamical histories. 
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